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Introductlon 

There  is  a  close  relationship  between  problems  of  structural 
optimization  and  the  analysis  of  media  with  microstructure.  The 
optimal  design  of  variable  thickness  plates  is  a  case  in  point:  for 
certain  problems,  plates  with  "stiffeners”  formed  by  rapid  thickness 
variation  can  be  stronger  per  unit  volume  than  any  traditional,  uniform 
or  slowly  varying  plates.  To  resolve  such  a  design  problem  one  must 
Introduce  a  "generalized  plate  model,"  representing  the  overall  effect 
of  a  microstructure  of  stiffeners  on  the  behavior  of  the  plate. 

One  idea  would  be  to  substitute  a  rapidly  varying  thickness 
function  into  the  fourth-order  equation  of  Kirchhoff  plate  theory  and 
perform  some  kind  of  "homogenization".  There  is,  however,  a  physically 
more  correct  approach:  it  appeals  directly  to  three-dimensional  linear 
elastostatics  on  thin,  rapidly-varying,  plate-like  domains.  There  are 
two  small  parameters  —  the  mean  thickness  e  and  the  length  scale  of 
thickness  variation  6  —  and  one  can  study  the  asymptotics  of  the 
solution  as  they  both  tend  to  zero.  This  was  the  focus  of  our  recent 
papers  [13,14],  We  showed  that  i£  makes  a  difference  which  parameter 
tends  to  zero  faster.  Use  of  the  Kirchhoff  plate  equation  with  a 
rapidly  varying  thickness  corresponds  to  the  case  e  «  6.  The  other 
extreme,  6  «  e,  corresponds  to  averaging  the  effect  of  the  thickness 
variation  first,  then  applying  Kirchoff  theory  to  the  resulting 
anisotropic  plate.  Intermediate  between  these  is  a  third  case,  e  -  6, 
which  has  no  such  simple  interpretation.  For  applications  to  optimal 
design  it  Is  natural  to  ask  which  alternative  gives  the  strongest 
structure,  and  that  was  the  focus  of  our  most  recent  paper  [15]. 

The  present  article  is  an  expository  review  of  this  work  and  its 
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relevance  to  optimization.  Special  attention  is  focused  on  plates  with 
"one  family  of  stiffeners,"  for  which  the  theory  is  relatively 
complete.  Much  remains  to  be  done  for  more  general  thickness 
variation;  various  open  questions  will  be  indicated  as  we  proceed,  and 
especially  in  section  6.  We  shall  refer  only  to  the  most  recent 
relevant  articles,  without  any  attempt  at  a  complete  survey  of  the 
extensive  literature.  More  references  on  homogenization  and  plate 
theory  can  be  found  in  [6,1 3. 2M]  and  an  extensive  bibliography  on 
structural  optimization  is  given  in  [33.  Recent  surveys  on  plate 
optimization  include  [2]  and  [20]. 


1 .  An  Optimal  Design  Problem 

Kirchhoff  plate  theory  models  the  behavior  of  symmetric,  variable¬ 
thickness  plates  under  transverse  loads.  It  specifies  the  vertical 
displacement  wQ  as  the  solution  of  an  elliptic  equation 


(1.1) 


32  32w0 

3xa3x0  3xY3xfi  ^  F 


on  the  midplane  domain  u,  with  appropriate  boundary  conditions  at  the 
plate  edges  3u.  The  tensor  MagYfi  relates  bending  moment  to  midplane 
curvature;  it  depends  on  the  plate's  thicknss  2h  and  on  the  constant 
elastic  moduli  of  the  material  from  which  the  plate  is  made, 
through  the  formula 


(1.2) 


Ma8Y« 


9 


where 


Ba$Y6  "  Bo$Y«  '  Bo633  BY633  /  B3333  * 


(The  Hooke's  law  tensor  is  assumed  to  satisfy  the  usual 

symmetries  -  Bjikt  “  Bki,ij»  2111(1  t0  have  the  mldPlane  aa  a  plane 

of  elastic  symmetry.)  For  an  isotropic  material,  B  is  given  by 


B1 1 1 1  m  B2222  m 


B1122  "  B221 1  m  Bv/(1~v^) 


B1212  *  B1 221  “  B21 1 2  *  B21 21  *  E/2(1+v)  , 


where  v  devotes  Poisson's  ratio  and  E  is  Young’s  modulus.  The  right 
side  of  (1.1)  is  the  load  per  unit  midplane  area. 

For  simplicity,  we  shall  discuss  only  plates  that  are  clamped  at 
the  edges;  this  means  that 


wQ  “  — —  •  0  at  3u  . 


The  principle  of  minimum  energy  gives  an  alternate  characterization  of 
Wq  as  the  minimizer  of 


(,-3)  S I  Mas«  5$: -  1F“ 

(i)  O  p  T  0  u 


Op 

in  the  Sobolev  space  Hc(w).  The  compliance  L  is  the  work  done  by  the 
load. 
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By  (1.3)*  it  has  the  variational  characterization 


(1.4) 
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For  a  given  load  F,  we  think  of  L  *  L(h)  as  a  functional  of  the 
thickness  h.  It  represents  an  overall  measure  of  the  plate's  rigidity 
under  F.  Therefore  it  is  natural  to  consider  the  problem  of 
optimization  for  minimum  compliance;  we  seek  to  minimize  L(h)  among 
all  plates  with  prescribed  volume  and  specified  minimum  and  maximum 
thickness,  l.e.  among  all  h  such  that 


(1.5)  h  e  L*((o)  ,  hmln  <  h  <  hfflax  and  j hdx  -  c  . 

(0 

It  Is  now  widely  recognized  that  for  some  choices  of  F  and 
hmax/hmin  thI3  optimal  design  problem  will  have  no  solution.  The 
difficulty  is  easy  to  understand  physically.  We  anticipate  that 
formation  of  "stiffeners"  by  means  of  an  oscillatory  thickness  could 
improve  the  strength  of  the  plate.  Since  tall,  thin  beams  are  stronger 
than  short,  fat  ones,  the  strength  should  increase  as  the  stiffener 
width  tends  to  zero.  If  there  is  no  optimum  scale  for  the  oscillation, 
then  there  will  be  no  optimal  h.  (A  more  precise  version  of  this 
argument  will  be  presented  in  section  2.) 

Numerical  manifestations  of  this  phenomenon  have  been  observed  in 
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[1,9].  Fop  certain  loads  F  and  sufficiently  large  ratios  hmax/brain» 
numerical  methods  for  minimizing  L(h)  are  seen  to  display 
instabilities.  The  computed  solutions  become  strongly  mesh-dependent, 
with  "stiffeners"  (oscillations  of  the  thickness  between  hmln  and  hmax) 
forming  on  the  same  scale  as  the  mesh  size. 

£ 

Mathematically,  the  point  is  that  L(h)  is  not  weak  lower 

semlcontinuous  on  the  space  (1.5)  of  admisible  h’s.  There  will  surely 

be  a  minimizing  sequence  { hn)  which  approaches  the  optimal  behavior, 

* 

and  (after  passage  to  a  subsequence)  it  will  have  a  weak  limit  h„,. 
But  he  compliance  can  jump  up  in  the  limit,  and  in  that  case  h*  will 
not  be  an  optimum. 

dearly  there  is  something  unsatisfactory  about  the  formulation  of 
a  design  problem  that  has  no  solution.  One  way  out  is  to  restrict  the 
design  space  by  imposing  a  pointwise  or  integral  bound  on  |Vh| 
(cf.  [5]).  The  other,  we  think  more  natural  alternative  i3  to  extend 
the  design  space  by  allowing  plates  with  stiffeners  or  rapidly  varying 
thickness  [^,10, 17].  This  entails  introduction  of  a  class  jB  of 
"generalized  plate-thicknesses"  and  an  extension  L  of  L  to  JJ  such  that 

( 1 ,6a)  For  each  fi  fc  the  generalized  compliance  L(fi) 
is  realizable  by  a  limit  of  ordinary  plates. 

In  other  words,  there  exists  a  sequence  {hnJ 
satisfying  (1.5)  for  which  C(B)  -  limn+00L(hn) . 

(1.6b)  The  functional  L  attains  its  minimum  value  on  £  . 

The  first  condition  assures  that  the  stiffeners  have  been  modelled 
correctly,  and  hence  that  the  underlying  problem  has  not  been  altered. 
In  particular,  it  implies  that  inf  L  -  inf  C.  The  second  condition 


-jVj'a'Vj  vVl  V 


says  that  the  class  /r  of  generalized  thickness  variations  is  "large 
enough".  It  promises  that  nothing  would  be  gained  (for  this  design 
problem)  by  considering  further  extensions  of  the  design  space. 

The  new  problem  of  minimizing  L  on  -cr  is  sometimes  called  a  full 
relaxation  of  the  original  design  problem.  (The  reader  is  warned, 
however,  that  this  term  is  used  slightly  differently  in  the  calculus  of 
variations,  for  example  in  [12].)  An  extension  to  some  intermediate 
class  of  plate  models  satisfying  (1.6a)  but  not  (1.6b)  could  be  called 
a  partial  relaxation.  Finding  a  partial  relaxation  requires  the 
correct  modelling  of  a  particular  class  of  plates  with  rapidly  varying 
thickness.  Finding  a  full  relaxation  is  more  difficult:  it  requires 
understanding  just  which  types  of  stiffeners  or  rapidly  varying 
thicknesses  can  occur  in  an  optimal  structure.  This  remains  in  general 
an  unsolved  problem,  but  the  easier  case  of  plates  with  a  "single 
family  of  stiffeners"  is  fairly  well  in  hand.  We  shall  discuss  it  in 
the  next  section. 

As  if  finding  a  relaxation  of  the  original  design  problem  were  not 
trouble  enough,  there  is  also  the  further  difficulty  of  its  relation  to 
three-dimensional  elasticity.  This  will  be  treated  in  sections  3  and 
4,  where  we  describe  a  class  of  three-dimensional  "plates"  with 
rap idly- varying  thickness  which  are  correctly  modelled  by 
homogenization  of  the  Kirchhoff  plate  equation  (1.1).  The  analysis 
shows,  however,  that  use  of  the  Kirchhoff  theory  above  represents  a 
loss  of  information:  plates  with  more  rapid  thickness  variation  require 
a  different  model.  Section  5  discusses  the  implications  of  this  for 
structural  optimization. 

Though  our  discussion  of  the  need  for  relaxation  has  focused  on 
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questions  of  existence,  the  relaxed  problem  is  as  important  for 
computation  as  it  is  for  the  theory.  Even  partial  relaxation  may  be 
advantageous  for  numerical  use.  Numerical  minimization  of  a  fully 
relaxed  C  will  be  free  of  the  instabilities  experienced  using  L;  also, 
experience  suggests  that  C  will  have  fewer  local  minima  than  L. 
Finally,  since  L  is  known  to  achieve  its  minimum,  one  can  obtain 
qualitative  information  about  extremal  designs  by  studying  the  first- 
order  optimality  conditions  for  L. 

2.  Rapid  Variation  and  Relaxation  of  the  Compliance  Functional 

In  order  to  relax  the  design  problem,  we  must  consider  how  rapid 
variations  in  h  affect  the  compliance.  There  is  a  general  theory  of 
homogenization  of  periodic  structures,  which  addresses  precisely  this 
sort  of  question  [6,24].  It  characterizes  the  vertical  midplane 
displacement  wQ  -  in  the  limit  as  the  length  scale  of  the  oscillation 
tends  to  zero  -  a  the  solution  of  (1.1)  with  a  new,  effective  rigidity 
fl  8yr .  The  limiting  compliance  is  correspondingly  / Fwq. 

The  simplest  case  is  that  of  a  plate  made  from  an  isotropic 
material  using  "stiffeners  in  the  X2  direction."  This  means  that  h  is  a 
function  of  x1  only,  independent  of  x2.  We  obtain  oscillations  on  a 
length  scale  6  by  taking  the  particular  form 

(2.1)  h5  -  H(x1  ,  x1  /6) 

where  H(x1 ,n1 )  is  periodic  in  the  second  variable  with  period  1 ,  and 
sufficiently  smooth  in  the  first  variable.  If  w^  is  the  solution  of 


1.1)  with  h  «  hg,  then  it  is  an  exercise  in  homogenization  to  see  that 
5  tends,  as  6  +  0,  to  the  solution  Wq  of 


1.2) 


a2  a2w0 

3xa3Xg  ^aBY6  axY3x6  ^  F 


Lth 


Hnn  ’  «<x,,-r3 

->  1  -V 

fl2222  “  I  E  +  |  ^4  H(Xl,.)"3  _1 

i-3)  _ 

R1 1 22  *  fl2211  “  |  ~~2  (H(x1»*)  3)  1 
5  1-v^ 

E  o 

H1 21 2  “  fl2112  ”  ^1 221  ”  fl2121  *  ffi+v)  H(x1’’)  * 

jre  H(x-,,.)3  denotes  the  average  of  the  periodic  function  H(x1f*)3 
Lth  respect  to  its  second  variable,  and  similarly  for  H(x.,,*)~3.  If 
-  H(x.|)  is  independent  of  »  1 . e.  if  there  is  no  rapid  variation, 
len  (2.3)  naturally  agrees  with  (1.2).  The  convergence  of  w^  towards 

Op 

j  is  in  the  weak  topology  on  H  ;  it  follows  that  the  compliances 
inverge 

! .4 )  L(hfi)  -  |  F  w6  - >  |  F  wQ  -  L  , 

u>  U) 

id  also  that  w^  -*•  Wq  uniformly  on  uj. 

To  see  the  advantage  of  rapid  thickness  variation,  we  consider 


jcillatory  perturbations  of  a  smoothly  varying  H: 
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olds  for  the  elastic  law  (5.5)  (and  more  generally  whenever 
1112  “  °)*  this  gives  an  example  in  which  Ma=1  -  Ma<1  is  indefinite, 
e  note  that  (5.5)  is  only  a  minor  modification  of  the  isotropic  law, 
hich  corresponds  to  the  choices 


5.8 ) 


Ev 


(1+v)(1-2v) 


2(1+v) 


We  hoped  at  first  that  this  would  lead  to  examples  of  practical 
ignificance.  If  the  difference  Ma“]1  -  Ma^]  1  were  large,  then  use  of 
ie  a  «  1  plate  (with  stiffeners  in  the  direction  of  greatest  bending) 
3uld  be  advantageous  for  some  design  problems.  However,  practical 
umerical  examples  of  (5.7)  are  exceedingly  hard  to  find.  For  the 
aometries  of  figures  1  and  2,  our  calculations  give  M3”1  <  Ma<1  for 
ie  elastic  law  (5.5).  The  best  example  we  found  -  and  it's  not  a  good 
ie  -  is  shown  in  Figure  3.  With  A  and  y  chosen  by  (5.8)  with  E  -  1.0 
id  v  -  0.25,  one  obtains  easily  that 


5.9) 
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The  a  <  1  model  does  not  depend  on  the  value  of  y'.)  On  the  other  hand 


5.10) 
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lere  £  is  the  strain  energy  of  the  cell  problem  (5.5)  -  which 
spends,  of  course,  on  y'.  For  the  isotropic  material  with  E  -  1.0  and 
■  0.25,  u’  is  equal  to  0.4;  for  the  calculations  reported  here  we 
sok  y'  *  5.0.  (Larger  values  of  y'  did  not  increase  the  significance 
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strictly  speaking  apply.  However,  the  arguments  presented  there  are 
easily  modified  to  include  this  case.) 

The  data  in  tables  1  and  2  naturally  satisfy  (5.2).  More 
interesting  is  the  observation  that  M3"1  and  Ma>^  are  quite  close, 
while  M32]2  is  much  greater  than  M32]2  in  each  case.  It  is  not 
surprising  that  figure  2  is  much  stiffer  than  figure  1 ;  we  understand 
that  N.  Olhoff  and  his  collaborators  are  currently  studying  the  use  of 
geometries  such  as  that  in  figure  2  for  compliance  optimization. 

Our  assertion  (5.4)  concerning  the  anisotropic  case  is  based  on  an 
explicit  counterexample:  for  an  elastic  law  of  the  form 

Biiii  -  ^  +  2m 

(5.5)  BilJJ  -  x  i  *  J 

b1212  "  m  »  Bl  31 3  *  b2323  “ 

with  p*  sufficiently  large,  and  for  a  thickness  profile  of  the  form 

(5.6)  H(n1 )  -  hQ(1  +  o*cos2iTrii  ) 

with  hQ  and  0  sufficiently  small,  we  showed  that 

(5.7)  M^],  <  Mf"],  and  m|^2  <  m|2J2  • 

Since  the  Inequality 


M 


a-1 

1212 


<  M 


a<1 

1212 
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(5.2)  Ma>1  <  M3"*1  <  Ma<^  for  an  Isotropic  elastic  law. 

The  left  Inequality  in  (5.2)  holds  even  for  anisotropic  materials: 

(5.3)  Ma>1  <  M3"1  in  general  (for  one  family  of  stiffeners). 

However,  the  right  inequality  in  (5.2)  can  fail  for  some  anisotropic 
laws  and  some  choices  of  the  profile  H: 

(5.4)  Ma>1  |  M3-1  in  general  (even  for  one  family  of  stiffeners). 

The  proofs,  presented  in  [15],  make  use  of  variational  principles  for 
each  of  the  three  effective  rigidity  tensors  Ma<1 ,  M3”1 ,  and  Ma>1 . 

More  quantitative  comparison  requires  numerical  calculation  of  the 
effective  rigidity  tensors.  Figure  1  shows  the  rescaled  cross-section 
(the  graph  of  H)  for  one  of  the  examples  presented  in  [13].  By  cutting 
along  the  horizontal  midline  and  glueing  together  opposite  ends  of  the 
stiffeners  we  get  a  new  cross-section  with  the  same  amount  of  material, 
Figure  2.  Tables  1  and  2  list  the  effective  rigidities  for  these 
geometries,  using  an  isotropic  material  with  Young's  modulus  E  -  1.0 
and  Poisson’s  ration  v  -  0.25.  The  cases  a  <  1  and  a  >  1  were  done 
using  explicit  formulas,  *#iich  are  easily  derived  as  in  [13]. 
Calculating  M3^  requires  finding  the  energies  of  a  pair  of  two- 
dimensional  cell  problems,  one  involving  plane  strain  and  the  other 
antiplane  shear.  This  was  done  using  the  FEARS  finite  element  code, 
developed  at  the  University  of  Maryland.  (Since  the  cross-section  in 
figure  2  is  not  of  the  form  (5.1),  the  analysis  of  [13.14,15]  does  not 


the  quadratic  forms  defined  by  their  rigidities  M  and  R.  We  say  that  M 
is  weaker  than  fl  if  M  <  fl,  in  other  words  if 


Maf5Y6  fca0  fcY 6  -  fla0Y<$  fcae  fcY 6 

for  every  symmetric  2x2  tensor  tfflg.  By  the  variational 
characterization  (1.H),  the  weaker  plate  has  the  greater  compliance 
under  any  load.  It  may  happen,  of  course,  that  neither  plate  is  weaker 
than  the  other;  in  that  case  the  ordering  of  the  compliances  will  be 
load-dependent. 

We  shall  suppose  for  simplicity  that  the  thickness  h  is  periodic 
and  depends  on  alone: 

(5.1)  h(x)  -  e  H(x1 /ea)  . 

As  remarked  in  section  2,  the  generalization  to  slowly  varying 
stiffener  geometry  or  direction  -  e. g.  to  axisymmetric  plates  with 
circumferential  stiffeners  -  is  Immediate.  Since  we  wish  to  compare 
the  different  scalings  a  <  1 ,  a  -  1 ,  and  a  >  1 ,  it  is  convenient  to  fix 
the  choice  of  H  in  (5.1).  The  corresponding  effective  rigidities  Ma<1 , 
M3”1 ,  and  Ma>1  are  as  discussed  in  section  notice  that  all  three 
plates  have  the  same  average  thickness  eH. 

Surprisingly,  the  ordering  of  the  strengths  of  these  structures 
depends  on  the  three-dimensional  elastic  material  used  to  make  them. 
For  plates  made  from  an  isotropic  elastic  material  using  one  family  of 
stiffeners,  the  a  >  1  scaling  is  weakest  and  the  a  <  1  scaling 
strongest: 
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regular  H  -  H(n),  taken  for  simplicity  not  to  depend  on  the  "slow” 
variable  x.  We  established,  among  other  things,  that  the  relative 
energy  error  is  of  order  /e  as  e  -*  0.  The  main  ingredients  of  the 
proof  are  a  pair  of  Integral  estimates,  an  averaging  lemma,  and  lots  of 
integration  by  parts.  The  first  integral  estimate  is  a  version  of 
Korn's  inequality  for  thin  domains  with  mean  thickness  e  and  thickness 
variation  on  length  scale  e,  making  explicit  the  dependence  of  the 
"constant"  on  e.  The  second  inequality  asserts  a  weak  form  of 
Kirchhoff's  hypothesis  for  the  elastostatic  displacement  ue,  as  a 
consequence  of  the  symmetries  of  the  problem  (3.1).  The  averaging 
lemmas  serve  to  replace  certain  rapidly-varying  expressions  by  their 
mean  values;  they  quantify  the  rate  at  which  a  periodic  function 
converges  (weakly)  to  its  mean  as  the  period  tends  to  zero.  A 
convergence  proof  could  probably  be  given  for  the  a  <  1  case  using 
similar  methods,  but  new  ideas  seem  required  to  handle  a  >  1.  A 
particular  stumbling  block  is  Korn's  inequality  for  Q£:  as  e  ♦  0,  the 
"constant"  blows  up  faster  when  a  >  1  than  it  does  for  a  <  1 . 

5.  Comparison  of  the  models  for  one  family  of  stiffeners. 

We  have  seen  that  three-dimensional  elasticity  supports  several 
different  models  for  thin  plates  with  rapidly  varying  thickness, 
depending  on  how  the  length  scale  of  thickness  variation  compares  with 
the  mean  thickness.  For  applications  to  structural  optimization,  it  is 
natural  to  ask:  which  scaling  gives  the  strongest  structure?  A  rather 
complete  answer  is  available  in  the  case  of  plates  with  "one  family  of 
stiffeners  in  a  specified  direction." 

One  way  of  comparing  the  "strength"  of  two  plates  is  to  compare 
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JL  Zij  (£a6)  -  0  in  Q(x) 

J 


(4.4)  - 


n3®iYot6vY  i"1 »2 


i-3 


on  3+Q(x)  , 


periodic  in  rj  . 


Here  Q(x)  denotes  the  rescaled  period  cell 
Q(x)  -  {g  :  | TI3 1  <  H(x;r))}  , 


Iij(£)  is  the  stress  EkS,(.£)  associated  to  the  strain 
Ek^(£)  -  +  and  v  13  the  outward  normal  to  Q;  notice 
that  x  enters  (4.4)  only  as  a  parameter.  The  formula  for  M3"1  is 

<"-5)  "Sels  ■  f  i*1  BaBf5  -  !  E14‘±°6>E1j<l1,S)'1''3  • 

-H 

where  as  usual  the  overbar  denotes  an  average  over  In  roost 
cases  the  actual  calculation  of  M3"1  must  be  done  numerically,  by 
solving  a  finite  element  or  finite  difference  approximation  of  (4.4). 

There  is  a  sense  in  which  the  a  -  1  model  includes  the  other  two. 
For  any  X  >  0,  one  can  apply  (4.4)  -  (4.5)  to  the  X-periodic  function 
Ha(x;3)  -  H(x;tj/X);  this  amounts  to  taking  6  -  Xe  in  (4.1).  The  tensor 
so  obtained  converges  to  Ma<1  (rigorously)  as  X  *  ®,  and  to  Ma>1 
(formally)  as  X  +  0;  see  [14]  for  details. 

The  convergence  of  three-dimensional  elastotatics  to  our 
generalized  plate  model  has  been  proved  with  mathematical  rigor  only 
for  the  Intermediate  case  a  -  1.  In  [14]  we  considered  a  sufficiently 
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Bapld  variation  of  the  thickness  of  the  "three-dimensional 
plate  is  correctly  modelled  by  homogenizing  the  Klrchhoff 
(A1 )  plate  equation  (only)  if  i_t  occurs  on  a  length  scale  larger 

than  the  average  thickness. 

There  Is  an  equally  simple  interpretation  for  Ma>1,  involving 
homogenization  of  a  rough  boundary.  Consider  the  system  of  elasticity 
on  the  domain  (4.1)  with  e  fixed,  as  6  ♦  0.  The  limiting  displacement 
will  solve  a  new  elasticity  problem  on  the  "smoothed"  plate  domain 

{(x,Xg):  x  e  «)  ,  |xg|  <  e  max^  H(x;tj)} 

with  non-constant,  effective  elastic  moduli  (see  the  appendix 

of  [13]  for  details).  Applying  Kirchhoff  plate  theory  to  this 
"smoothed"  structure  corresponds  to  taking  the  limit  e  -*•  0,  and  it 
yields  precisely  the  rigidity  tensor  Ma>1.  (This  calculation  remains 
at  present  merely  formal.  The  homogenization  of  rough  boundaries  was 
made  rigorous  for  a  scalar  equation  in  [7],  but  the  system  of 
elasticity  presents  additional  difficulties.) 

The  intermediate  case  a  -  1  has,  unfortunately,  no  such  simple 
interpretation.  The  tensor  M3"1  depends  on  auxiliary  functions 
i°*(  x;rj),  obtained  by  solving  the  elastotatlc  boundary  value  problem 


-16- 


We  are  convinced  that  (4.1)  and  (4.2)  represent  no  less  of  generality 
toward  the  problem  of  characterizing  the  limiting  compliances. 

The  asymptotics  of  ue  as  e  •*  0  can  be  studied  using  the  method  of 
multiple  scales  [13.14].  (A  closely  related  problem,  involving 
smoothly  varying  thickness  but  rapidly  varying  elastic  moduli,  was 
analyzed  simultaneously  and  independently  by  Caillerie  [8].)  The  main 
conclusion  is  that  the  energy  of  the  elastic  displacement  ue  approaches 
that  of  a  fourth-order  equation 


(4.3) 
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Since  the  compliance  is  -2  times  the  energy,  the  compliances  also 
converge.  The  solution  w  of  (4.3)  is  the  limiting  vertical 
displacement;  therefore  MagY$  represents  an  effective  rigidity  tensor , 
relating  midplane  curvature  to  bending  moments  in  the  limit  as  e-*0. 
It  depends  not  only  on  H(x;rj)  -  corresponding  to  the  geometry  of  the 
variation  -  but  also  on  the  choice  of  scaling,  specifically  on  whether 
a<1,a-1,ora>1.  Formulae  are  given  in  [13]  for  each  of  the 
rigidities  Ma<1 ,  M3"1,  and  Ma>1,  in  terms  of  the  solutions  of  certain 
periodic  boundary-value  problems. 

a/1 

The  tensor  MaNI  has  a  simple  interpretation:  it  is  precisely  what 
one  obtains  by  homogenizing  the  plate  equation  (1.1).  Thus  for  "one 
family  of  stiffeners"  and  an  isotropic  elastic  law,  Ma<1  is  given  by 
(2.3).  In  particular,  this  analysis  yields  an  answer  to  our  first 
question  (Q1): 


EBBH 
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question  is  that  (2.2)  correctly  models  plates  whose  thickness  varies 
on  a  length  scale  that  is  large  compared  with  the  mean  thickness  but 
small  compared  with  the  plate  diameter.  We  do  not  have  a  complete 
answer  to  (Q2);  however,  we  are  able  to  compute  the  effect  of  thickness 
variation  on  a  length  scale  comparable  to  or  shorter  than  the  mean 
thickness.  These  lead  to  different  results  than  those  predicted  by  the 
Kirchhoff  theory. 

4.  Summary  of  our  generalized  plate  models. 

The  first  important  step  towards  answering  the  questions  Q1  and  Q2 
is  to  understand  the  influence  of  rapid  thickness  variation  on  the 
(rescaled)  compliance  of  a  "three-dimensional  plate", 

I  |  F | n| | u|  -  e~3(e3  |  F | n| | u|  -  j  o(ue) *e(ue ) ]  , 

3+Qg  d+Qg  Og 

where  3+Q£  denotes  the  upper  and  lower  surfaces  of  the  plate.  We  shall 
assume  that  the  thickness  is  locally  periodic  with  period  6  -  6(e);  in 
other  words 

(4.1 )  -  {(x,x3):  x  e  «,  |x3|  <  eH(x;x/6)}  , 

where  H(x;g)  is  periodic  in  its  second  variable  with  period  1.  Though 
the  dependence  of  6  on  e  could  in  principle  be  arbitrary,  it  is 
convenient  to  choose 


(4.2) 
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0  <  a  < 


(3.2) 
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where  wQ  solves  (1.1). 

The  convergence  (3.2)  holds  provided  that  h  and  F  are  fixed  and 
sufficiently  smooth.  However,  we  have  seen  that  relaxing  the  design 
problem  requires  introduction  of  "generalized  plate  models"  based  on 
rapidly  varying  thicknesses.  Since  this  appears  to  violate  the 
hypotheses  on  which  (3.2)  is  based,  it  is  important  to  ask: 

What  is  the  relation  between  the  homogenized  plate 
(Q1 )  equation  (2.2),  its  relaxed  compliance  E,  and 

solutions  of  the  equations  of  three-dimensional  elasticity? 

This  leads  naturally  to  a  second,  much  harder  question: 

Consider  linear  elasticity  on  the  three  dimensional  plate 
domain  |x3|  <  e  he(x1tx2)»  (x1tx2)  e  to,  with 
vertical  load  e^F  per  unit  midplane  area.  Can  one 
(Q2)  characterize  ( in  a  simple  way)  the  lowest  limiting 

point  that  may  be  obtained  for  the  compliances,  as 
e  ->0,  by  letting  he  vary  with  e  subject  to 

^min  1  -  hmax*  jhe  “  c? 

u 

These  issues  are  the  focus  of  our  recent  work  [13,14,15]  which 
will  be  summarized  in  the  next  two  sections.  Our  answer  to  the  first 


also  [233.  One  route  to  a  full  relaxation  involves  seeking  "optimal 
bounds"  for  the  effective  rigidity  of  a  plate  in  terras  of  hp,^,  h,,,^, 

and  the  mean  thickness  h;  another  involves  the  study  of  lower 

semi continuity  and  the  calculus  of  variations.  Both  approaches  have 
successfully  treated  related  problems  [12,17,19,21]  and  it  seems 
reasonable  to  hope  for  further  progress  soon. 

3.  An  Apparent  Physical  Contradiction. 

The  primary  justification  for  the  plate  equation  (1.1)  is  that  it 
follows  from  the  equations  of  linear  elasticity  on  the  three- 
dimensional  plate  domain  in  the  limit  as  the  thickness  tends  to  zero. 

Specifically,  let  Q£  denote  the  plate  domain 

Pe  -  {(x1,x2,x3):  (x1tx2)  e  u  ,  |xg|  <  e  hCx^Xg)}  , 

and  let  ue  be  the  linearly  elastic  displacement,  satisfying: 
div  o(ue)  -  0  in  flE 

(3.1)  o(ue)*ne  -  1  e3(oto,F|n|| )  for  |xg|  -  +  eh 

ue  -  0  for  ( x1 ,  x2 )  e  3w  . 

Here  o(ue)  denotes  the  stress  ek£  associated  to  the  strain  ek4  - 

i(3uk/3x^  ♦  3u|/3xk),  and  ne  is  the  outward  unit  normal  to  fle.  Using 
either  the  dual  variational  principles  [18]  or  a  direct  asymptotic 
expansion  of  ue  [11],  one  can  show  that  the  (rescaled)  energy  converges 
as  t  *  0  to  that  given  by  the  Kirchhoff  model: 


following  the  ideas  of  CIO],  The  spaced  of  "generalized  thickness 
variations"  consists  of  pairs  R  =  (g,0).  The  first  element  is 
restricted  by  h,^  <  g(x.,)  <  hmax;  it  represents  the  "minimum  height  of 
the  stiffeners."  The  second  is  constrained  by  0  £  e(x.|)  <  1,  and  plays 
the  role  of  the  stiffener  density.  Both  g  and  6  are  assumed  to  be 
measurable,  but  not  necessarily  continuous.  The  rigidity  tensor  W0gY6 
corresponding  to  fi  is  defined  by  (2.3)  with 


H(x1  Jni)  - 


max 


g(Xl  ) 


o  <  n1  <  e(x1 ) 
e(x.,)  <  ni  <  i  » 


and  the  corresponding  compliance  L(R)  is  obtained  by  solving  (2.2). 

This  defines  at  least  a  partial  relaxation,  by  virtue  of  (2.4).  We 
believe  that  it  is  the  full  relaxation,  though  this  has  not  yet  been 
established  with  mathematical  rigor.  Optimality  conditions  for  the 
analogous  relaxation  of  axisymmetric  plates  are  presented  in  [10];  they 
determine  the  circumstances  under  which  "stiffeners”  occur  in  an 
optimal  structure. 

The  situation  becomes  much  more  complicated  if  more  general 
thickness  variation  is  allowed.  For  "two  or  more  families  of 
stiffeners"  the  effective  rigidity  can  only  be  expressed  in  terms  of 
the  solutions  of  certain  fourth-order  equations  on  the  period  cell  (see 
e.  g.  [13]).  For  a  "single  family  of  stiffeners"  w'th  unspecified 
direction,  there  is  the  possibility  of  the  direction  itself  becoming 
oscillatory.  And  there  are  more  general  geometries  to  consider, 
analogous  to  he  composites  of  rank  two  or  more  considered  in  [16],  see 
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I 


n  -  -j 


32w, 


^o6Y6  3x  3xe  3xY3xfi  +  2  i  FwA 

u  Ot  p  T  0  Oj 


f  32wx  *2"\  f 

<  -  M  x — r—  r — 7-—  ♦  2  Fw, 

J  o6Y6  3w  3x0  dXydXr  J  * 

(i)  Ot  p  T  0  (t) 

<  L(H)  , 


and  therefore  the  oscillations  have  improved  the  compliance. 
Intuitively,  deformations  of  the  desired  type  are  expected  (at  least 
away  from  the  edge  of  the  plate)  whenever  F  oscillates  rapidly  enough 
with  respect  to  x2. 

The  preceding  discussion  is  easily  extended  to  a  "single  family  of 
stiffeners"  with  smoothly  varying  profile  and  direction.  This 
corresponds  to  choosing 


h6  -  H(x;  e(x)*x/6)  ,  |e(x)|2  -  1  , 


in  place  of  (2.1).  The  stiffeners  are  then  orthogonal  to  the  field  of 
unit  vectors  e(x);  their  profile  is  determined  by  H(x1tx2;n)»  which 
should  be  periodic  in  the  real  variable  n  with  period  1 .  Since 
homogenization  is  local,  the  effective  rigidity  at  any  x  e  u>  will  be 
given  by  (2.3)  in  the  orthogonal  coordinate  system  which  takes  the  x^ 
axis  parallel  to  e(x).  The  case  of  axially  symmetric  plates  with 
circumferential  stiffeners  is  an  especially  natural  one.  It  was 
treated  in  [9,10],  and  rapid  thickness  variation  was  found  to  be 
advantageous  for  loads  of  the  form  F  -  cos  ke  when  k  *1. 

Returning  for  simplicity  to  plates  with  stiffeners  in  the  x2 
direction  (l.e.  h  a  function  of  x1  alone),  we  formulate  a  relaxation  by 
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H^(x.j , Tii  ^  ™  H(Xi )  +  A^(tii  )  , 


[♦(n1 )dn1  -  0 


If  <  H  <  hmax  then  the  same  will  be  true  for  when  X  is  small 
enough.  Calculation  gives  that 


—  H?-  0  ,  —  H3  ■  6  I 

«  1  d»2 


£  (H-3,-1  -  0  ,  ±-  -  -12  H 


at  X  *  0.  We  see  from  (2.3)  that  a  small,  oscillatory  perturbation 
always  decreases  and  whlle  *t  increases  R1212.  In  the 
physical  range  0  <  v  <  1,  0<E,  it  also  increases  ^222*  7,113  la 
entirely  reasonable,  since  stiffeners  should  resist  twisting  and 
lengthwise  bending,  but  should  be  rather  weak  under  bending  in  the 
orthogonal  direction. 

We  assert  that  for  suitable  loads  F,  the  compliance  is  decreased 
by  such  perturbations.  Indeed,  let  be  the  effective  rigidity  of 
the  (oscillatory)  perturbed  geometry,  that  of  the  (smooth)  H,  and 
suppose  that  the  solution  w^  of  (2.2)  satisfies 


I2  «  I  -V  ♦  I 


*  12 


3x1 3x2 


Then  the  perturbed  compliance  C  satisfies 


of  the  findings,  but  only  served  to  increase  the  numerical  error.)  The 
value  of  ^  computed  by  FEARS  was  0.7119  *  10”\  using  574  elements  in 
a  subdivision  based  on  Just  1/4  of  the  unit  cell.  The  code  has  a 
built-in  error  estimator,  which  says  in  this  case  that  the  energy  is 
off  by  at  most  0.95$.  Extensive  practical  experience  with  FEARS 
indicates  that  the  true  error  exceeds  the  estimated  error  by  at  most 
75$.  We  are  thus  convinced  that  the  correct  value  of  in  (5.10)  is 
bounded  above  by  .7244x10”^,  which  yields 

Mfm  >  0.192951  *  10-2  . 

In  view  of  (5.9),  this  gives  a  numerical  example  of  (5.7)  -  but  hardly 
one  of  any  practical  significance! 

The  implications  of  these  results  for  structural  optimization  are 
clear:  for  the  design  of  plates  made  from  an  isotropic  material  using 
one  family  of  stiffeners  with  a  specified  direction,  attention  may  be 
restricted  to  the  a  <  1  model  (obtained  by  homogenizing  the  Kirchhoff 
plate  equation).  For  plates  made  from  an  anisotropic  material  this  is 
not  true,  but  our  numerical  experimentation  suggests  that  even  so 
little  will  be  lost  in  practice  by  a  restriction  to  Kirchhoff  theory. 

Care  is  advised  in  extending  these  conclusions  to  more  general 
situations,  such  as  plates  with  two  or  more  families  of  stiffeners.  We 
proved  in  [153  that  Ma<1  is  strongest  and  Ma>1  weakest  if  the  elastic 
law  satisfies  Bag33  *  0  -  which  includes  the  isotropic  law  with 
Poisson’s  ration  v  set  equal  to  zero.  But  in  general  the  relative 
strengths  appear  to  depend  on  both  the  elastic  moduli  and  the  form  of 
the  thickness  variation. 


6.  Directions  for  the  future 


We  indicate  some  of  the  areas  in  which  further  work  is  needed. 

a)  The  correctness  of  our  a  >  1  model  as  a  limit  of  three- 
dimensional  elasticity  has  yet  to  be  proved. 

b)  We  still  hope  for  an  example  (perhaps  using  two  families  of 
stiffeners)  where  the  a  -  1  model  leads  to  a  significantly  stronger 
structure  than  does  homogenization  of  the  Kirchhoff  theory. 

c)  Section  2  described  a  relaxation  of  the  compliance 
optimization  problem  for  Kirchhoff  theory  with  "one  family  of 
stiffeners  in  a  specified  direction."  It  remains  to  prove  an  existence 
theorem,  l.e.  to  establish  that  this  is  the  full  relaxation.  The 
corresponding  problem  without  the  restriction  to  "one  family  of 
stiffeners"  is  more  difficult.  A  relaxation  (partial  or  full?)  has 
been  proposed  in  [23];  we  understand  that  Lurie  and  Cherkaev  also  have 
made  progress  in  this  direction. 

d)  We  have  shown  that  Kirchhoff  theory  suffices  for  the 
optimization  of  plates  made  from  an  isotropic  material  with  one  family 
of  stiffeners.  This  result  seems  likely  to  fail,  however,  for  more 
general  geometries.  Therefore  the  best  relaxation  would  be  one  based 
not  on  Kirchhoff  theory  but  instead  on  threecdimensional 
elasticity  -  l.e.  on  our  generalized  plated  models. 

e)  We  have  discussed  only  the  compliance  optimization  of  linearly 
elastic  plates.  In  plasticity,  the  analogous  problem  is  to  maximize 
the  limit  multiplier  of  a  given  load.  Rapid  thickness  variation  arises 
naturally  in  that  context,  too,  and  relaxed  formulations  have  been 
proposed  in  [22,25].  It  remains,  however,  to  prove  that  the  models 
used  there  correctly  represent  the  asymptotic  behavior  of  three- 


dimensional  plastic  structures.  It  also  remains  to  check  whether  a 
full  relaxation  has  been  achieved. 
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a<1 

a-1 

a>1 

M1 1 1 1 

.015 

.012 

.011 

M1122 

.004 

.003 

.003 

M2222 

• 

U> 

CO 

4r 

.334 

.334 

M1 212 

.133 

.006 

.004 

Table  1 :  Effective  rigidities  for  figure  1 . 


a<1 

a-1 

a>1 

M1 1 1 1 

.777 

.687 

.678 

M1 122 

.194 

.172 

.169 

M2222 

.851 

.845 

.844 

M1 21 2 

.321 

.262 

.254 

Table  2:  Effective  rigidities  for  figure  2. 
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